Quantum interference in timed Dicke basis and its effect on bipartite entanglement 
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We analyze the effect of position dependent excitation phase on the properties of entanglement between two 
qubits formed in atomic systems. We show that the excitation phase induces a vacuum mediated quantum inter- 
ference in the system that affects the dynamical behavior of entanglement between the qubits. It is also found 
that the quantum interference leads to a coherent population transfer between the symmetric and antisymmetric 
states which can considerably modify the dynamics of two-qubit entanglement and can even prevent finite time 
disentanglement (sudden death) under certain conditions. 
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I. INTRODUCTION 



Quantum interference (QI), an intriguing consequences of 
the superposition principle has led to numerous fascinating 
phenomena such as coherent population trapping JTt], las- 
ing without inversion electromagnetically induced trans- 
parency H, and quantum entanglement 04al- The applica- 
tion of QI in generation of bipartite entanglement both in dis- 
crete 10, H] and continuous-variable JUGS] settings has been 
the focus of current investigation. Note that bipartite entan- 
glement involving two atoms, extensively used for implemen- 
tations of various quantum information protocols lUlUHjll. i s 
known to be quite fragile in the face of decoherence ll 1 ll 1 nil . 
In view of this, in the past few years considerable effort has 
been devoted to the study of dynamical aspect of two atom en- 
tanglement in presence of decoherence iflia - ual . In one such 
study [ 18J it was found that in contrary to the adverse effect 
of spontaneous emission on atomic entanglement |24j], coop- 
erative spontaneous emission in two atom systems can gen- 
erate entanglement among the atoms. It is worth mention- 
ing here that the problem of cooperative spontaneous emis- 
sion first addressed by Dicke 12711 is known to exhibit sev- 
eral counter-intuitive phenomena in two atom systems 112811 
namely, directed spontaneous emission |29ll. Lamb shift IB0I1 . 
single photon Dicke superradiance 113111 and others. In re- 
cent times, with the discovery of atom like behavior of semi- 
conductor quantum dots 11331 — [351] an d their utilization towards 
solid state quantum computing J33l - [36ll . we have a new class 
of systems where the phenomenon of cooperative spontaneous 
emission can be of immense importance from the context of 
quantum information sciences. 



Recently, Ooi et al. 113711 studied the effect of position 
dependant excitation phase on the population dynamics, in- 
tensity and spatial and angular correlations for two two-level 
atoms interacting via their dipoles. The results show that the 
excitation phase considerably modifies the dynamics of the 
system. Later, Das et al. 13811 investigated the effect of the 
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position dependent excitation phase on the Dicke cooperative 
emission spectrum. A strong quantum correlation among the 
atoms was reported in presence of the excitation phase. This 
was attributed to a vacuum mediated QI generated in the two 
atom system in presence of the position dependent excitation 
phase. The result of [38] qualitatively indicates that the spa- 
tial variation of the excitation phase can affect the generation 
and evolution of entanglement in the system. It may be added 
that, such vacuum mediated QI and its effect has been ear- 
lier studied in atomic systems lEM I39l - l4ltl . Further, a recent 
work has predicted how one can use such QI to protect bipar- 
tite entanglement l42tl . While these earlier works utilize the 
quantum interference that comes about due to the configura- 
tion of the atomic system, we are motivated at studying the 
effect of quantum interference induced by the position depen- 
dent excitation phase. 

To understand the effects of such QI on the two atom en- 
tanglement, we in this paper perform a systematic study of 
the time evolution of entanglement measure for two strongly 
dipole coupled atoms undergoing a cooperative spontaneous 
emission. We consider various initial quantum states in which 
the two atoms can be prepared and explore the effects on the 
dynamical behavior of entanglement that results as a conse- 
quence of the quantum interference. We explicitly take into 
account the position dependent excitation of the atoms by in- 
troducing timed Dicke basis ll29ll . It is important to understand 
that the entanglement in a two atom system crucially depends 
on the cooperative decay rates, the initial conditions, and the 
dipole-dipole interactions [21], all of which gets modified due 
to the quantum interference. It is worth mentioning here that 
QI arising from position dependent phase in such timed Dicke 
basis was explored in a recent study in context to population 
dynamics and photon correlation studies in two atom systems 
n37ll . We however in this current work are interested in investi- 
gating the effect of such QI on the entanglement of two atoms. 
For instance, for the system initially prepared in the symmet- 
ric timed Dicke state, a coherence between the symmetric and 
antisymmetric states is dynamically generated as a result of 
the QI between the two pathways to the ground state. This 
coherence leads to considerably slow decay of entanglement. 

The organization of the paper is as follows. In Sec. II we 
discuss our model and write down the dynamical equation for 
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our system using a master equation approach. Then in Sec. Ill 
we discuss the entanglement measure and derive generalized 
analytical expressions in the timed Dicke basis for the two 
atom system. In Sec. IV, we then consider two initial condi- 
tions: for the atoms prepared in pure states and mixed state 
and show explicitly that the vacuum mediated QI induced by 
excitation phase can lead to considerable modification of two 
atom entanglement behavior. We provide analytical and nu- 
merical results in support of our propositions. Finally, we 
summarize our results in Sec. V. 



II. MODEL AND EQUATIONS OF EVOLUTION 

We consider a system of two qubits formed by the excited 
states |ej) and ground states \gj) (i = 1,2) of two identical 
two level atoms. The qubits are fixed at positions n and r 2 
and the inter-atomic distance is less than the wavelength of 
the radiation field, A. We further assume that the qubits are 
coupled to one another by a dipole-dipole interaction and are 
coupled to the environment via an interaction with a common 
vacuum reservoir. The time evolution of the density operator 
for such a two-qubit system can be treated in a master equation 
framework and is given by 12811 
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where uiq is the atomic transition frequency, erf = (ajai — 
<7i<j\)/2 is the energy operator with <r](oi) being the raising 
(lowering) operator for ith atom, fly and 7^ for i 7^ j are 
respectively, the dipole-dipole interaction term and the coop- 
erative decay rate given by 
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where 2j = 2jn — 2j 2 2 — 2\p eg \ lu 3 /Sirsohc 3 is the spon- 
taneous decay rate of the individual qubits. p eg is the dipole 
moment, ko — 2tt/X with A being the wavelength of the emit- 
ted radiation and 9 is the angle between the direction of the 
dipole moment and the line joining the ith and the j th qubits, 
and rij = \t{ — Tj \ is the interqubit distance. In this paper we 
assume that the orientation of the dipole moment is random 
and hence Eqs. (O and ([3]) simplifies considerably and take 
the form 



fii, = -7cos(fc rjj)//c 7\,-, 
jij = 7sin(fe ry)/feory. 



(4) 
(5) 



We next consider the preparation of initial state of the 
qubits. For this purpose we assume that the qubits interacts 
with a very weak laser field (almost at a single photon level) 
propagating with a wave- vector ko. The interaction with the 
weak field can lead to a resonant single photon absorption pro- 
cess. It is important to note that we consider the direction of 
the wave-vector to be different to that of the inter-qubit axis. 
This thus generate a position dependent excitation phase of 
the qubits when ever a photon is absorbed. The excitation pro- 
cess, with the laser field treated classically and in the rotating 
wave approximation, can be described by the Hamiltonian 



V = -m^frfe* ** e- i ( v °- u °) t + adj.), 



(6) 



where SI = d ei9l - £/h = d e2ff2 -£/His the Rabi frequency and 
vq is the angular frequency of the incident radiation. Note that 
the position dependent phase factors in the Hamiltonian would 
substantially affect the dynamical behavior of the correlation 
in the two qubit system. We except that this in turn will lead 
to modification of entanglement among the qubits. The inves- 
tigation of any such modification in the entanglement feature 
is the key focus of this paper. In order to investigate the ef- 
fect of position dependent excitation phase on the dynamics it 
proves to be convenient to work in a basis defined by the phase 
factors. Such a basis was introduced in Ref. [29] in context 
to directed spontaneous emission from an ensemble of atoms 
and is also known as the timed Dicke basis. To this end, for 
our system of two qubits there are four timed Dicke states: 
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In terms of this basis the equations of evolution for the ele- 
ments of the density operator are: 
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FIG. 1: (Color Online)Energy level diagram for two two-level atoms in bare basis (a) and in the timed Dicke basis (b). The frequency shift 
A = f2i2 cos <p occurs as a result of dipole-dipole interaction between the two atoms. The collective states js) and \a) decays at a rate of T+ 
and r_ respectively, where r± = 2(7 ± 712 cos ip). 



p ss = - 2(7 + 712 cos ip)p ss - i sin (^(712 + iO, 12 )pas 

+ iSin i/3(7l2 -iVti 2 )p sa +2(7 + 712 COS</?)p ee , 

(lie) 

Paa = - 2(7 - 712 COS (f)p aa ~ 1 Sin tp(~fl2 ~ i^l2)Pas 

+ i sin (,9(712 + iVL 12 )p sa + 2(7 - 7i 2 cos^)p ee , 

(Hf) 

pas = - 2(7 - i£li2 cos(^)p as + i sin 1^(712 + £fii2)p ss 

+ i sin (^(712 - i£li2)Paa - simppee, (llg) 

Pgs = - [7 + 712 COS (£> - i{u Q + n i2 COS ^)]P flS 
+ i sin (^(712 - i^l2)Pga 

+ 2(7 + 7i 2 cos(^)p se + 2i7i 2 sin^p ae , (llh) 

Pga = - [7 - 712 COS <p-i(tJ - ^12 COS f)])p ga 

- i sin (^(712 - ifiia)p 9S 

- 2(7 - 712 cos <p)p ae + 21712 sin (pp se , (Hi) 

/9 5g =2(7 + 712 COS </2)p ss + 2(7 - 712 COS tp)p aa 
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where y> = ko • (rj — ) = cos f with ^ being the angle 
between the laser propagation direction and the line joining 
the two atoms. 

Inspection of Eqs. ( II Ibi and ( 111 cb shows that the presence 
of atom-atom interaction gives rise to collective level shift The 
level shift arises due to the atom-atom interaction only occurs 
in states \s) and \a). The other collective states |e) and \g) do 
not see any level shift due to this interaction as per Eq. (lid). 
That is, the energy difference between the |e) and \g) remain 
2u! . The state \s) is shifted up while the state \a) is shifted 
down by an equal amount A = fii2 cos <p from the single 
photon resonance line as shown in Fig. [lj>. It is interesting to 



note that one can manipulate the level shift by only orienting 
the laser field appropriately with respect to the line joining the 
two atoms. For example, ip — tt/2, i.e., when the angle be- 
tween the laser propagation direction and the line joining the 
two atoms is £ = tt/3 and the interatomic distance equal to 
half of the radiation wavelength, ri2 = A/2, the level shift 
vanishes. Thus it is possible to control the level shift by ap- 
plying a laser field in a particular direction without turning off 
the dipole-dipole interaction. 

Further, we note that the transition probability from the ex- 
cited state |e) to the one photon states, \s) and \a), is the sum 
of the probability of each transition. Since it is the probabil- 
ity, and not the probability amplitudes that adds up we don't 
expect quantum interference phenomenon to occur. However, 
the transition probability from the one photon states, \s) and 
I a) to the ground state \g) is obtained by squaring the sum of 
the amplitude of each transition. When there is a coherence 
between the two states (|s) and |a)), this can lead to quantum 
interference yielding coherent population transfer between | s) 
and I a). Indeed, the populations in \s) and \a) is coupled to 
the coherence p as as per Eqs. dl le| >-( p~Tg] >. It is worth to note 
that this coupling disappears when the direction of propaga- 
tion of the laser field is perpendicular to the interqubit axis 
£ = 7r/2 {ip = 0). Therefore, we see that in the presence of 
a position dependent excitation phase <p quantum interference 
is induced in the system. In this paper we hence explore to 
what extent the quantum interference developed in the system 
affects the dynamical properties of the bipartite entanglement 
between the two qubits. 



III. ENTANGLEMENT MEASURE 

In general a state of a quantum system is said to be entan- 
gled when the density operator of the composite system can- 
not factorize into that of the individual subsystems. There are 
several entanglement measures for two-qubit system in the lit- 
erature. However, we use the concurrence, a widely used en- 
tanglement monotone, for our purpose. The concurrence, first 
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introduced by Wooters B43H . is defined by 



C(t) = max(0, V A i - V A 2 



A 4 ), (12) 



where Ai > X 2 > A3 > A4. are the eigenvalues of the 
matrix pp in which p — a y ® <j y p*a y ® a y with a y being the 
Pauli matrix. The concurrence takes values ranging from to 
1. For maximally entangled state C(t) = 1 and for separable 
state C(t) = 0. 

In general, for a dissipative system, without any external 
driving field, the density matrix of the qubits system has the 
form 
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in the following basis set 
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Note that for a quantum state initially prepared in a block form 
of ( fT3l . the time-evolved density matrix will have the same 
block form, i.e., the zeros remain zero and the nonzero com- 
ponents evolve in time ll2ll [26| • We next proceed to calculate 
the concurrence for the qubits system initially prepared in the 
form of ( fT3] l. To do so, one has to determine the matrix p 
in the basis where p is expressed. Using the definition of the 
density matrix p, we obtain 



Pit) 




(15) 



Thus the square root of the eigenvalues of the matrix pp are 

= {VP22P33 ± |P23|, \/PllP44 ± \Pli\}- 



(16) 



There are two possible expressions for the concurrence, de- 
pending on the values of the eigenvalues. The first case is that 
when |yc>23 1 + y 1 P22P33 be the largest eigenvalue. This leads to 
a concurrence 



Ci{t) = 2(1/3231 - y/pllm)- 



(17) 



While if 1^14 1 + \/PiiPii is the largest eigenvalue then the 
concurrence takes the form 



C 2 (t) = 2(\p U \ ~ VP22P33)- 



(18) 



Depending on the initial condition used, one of the con- 
currence expressions suffices to measure the entanglement 
present in the qubits system. Further, inspection of (TTTb and 
(fT8l shows that C\ (t) would be positive and hence the mea- 
sure of entanglement when the one photon coherence is larger 



the square root of the product of the populations in the excited 
and ground states. On the other hand, for C 2 (t) to be a mea- 
sure of entanglement for the system the two photon coherence 
should be greater than the square root of the product of the 
population in one photon excited states. 

In order to gain insight into the physics it is convenient to 
express the concurrences in terms of timed Dicke basis intro- 
duced earlier. To do so, one has to apply a unitary transforma- 
tion UpW on the density matrix given by (fT3l . The matrix U 
is given by 



U 



( 1 \ 

' J x - 

u V2 V2 

— — 

U V2 V2 

V 1 / 







(19) 



The elements of the density matrix UpU^ is related to that of 
pby 



(P22 + P33 - (P23 + P32)) 
(P22 + P33 + P23 + P32) 
33 + P23 ~ P32) 
33 ~ (P23 - P32))- (20) 



Therefore the concurrence can be expressed in terms of the 
timed Dicke basis as 
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= Pll 


Peg 


= Pli 






Paa 
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C(t)=max(0,Ci(*),C 2 (f)), 



(21) 



where 



C\{t) = V (Pss - Paa) 2 + 4[Im( ( o QS )] 2 - 2^/p ee p gg (22) 



C 2 {t) = 2\p eg \ - ^f(p ss +p aa ) 2 +4[Re(p as )] 5 



(23) 



This expression for concurrence will be used in the following 
section to study the dynamical evolution of entanglement in 
the two-qubit system by considering various initial conditions. 



IV. ENTANGLEMENT DYNAMICS OF TWO IDENTICAL 
QUBITS 

Using the entanglement measure introduced in the previ- 
ous section we study the effect of the position dependent ex- 
citation phase by considering pure and mixed state as initial 
conditions. 



A. Initial pure states 

In the two-qubit system one might consider a pure sepa- 
rable or entangled state as an initial condition. For instance, 
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FIG. 2: (Color online) A scheme illustrating a proposed method to prepare symmetric timed Dicke state |s). A similar scheme has been 
proposed to excite one atom in a cloud of N atoms |29]. The nonlinear crystal (NLC) down converts the a pump photon into signal-idler pair. 
A click on detector (Dl) indicates generation of the pair and hence no click on the second detector (D2)-assuming a perfect detector-means 
the other photon (ko) conditionally excite one of the atoms. 



for pure separable state, one can take the two atom excited 
state, |e). Even though this is unentangled state at the ini- 
tial time, the interaction of the atoms with the environment 
leads to weak transient entanglement 11181 12111 . The effect of 
quantum interference induced by position dependent excita- 
tion phase is unimportant in this case and thus we rather focus 
on pure entangled state as an initial condition. 

We take the initial state of the two-qubit system to be the 
symmetric state \s). This state is a pure maximally entan- 
gled state and can be prepared using correlated pair of pho- 
tons generated from a parametric down conversion process in 
which one of the photons is sent to a detector (Dl) and the 
other is directed towards the atoms. A click on the detector 
(Dl) tells us that the other photon is sent to the atoms. If the 
second detector (D2) registrars a count then no atom is ex- 
cited. However, if Dl shows a click and D2 does not then 
we know that one of the atoms is excited, but we don't know 
which one (see Fig. [2}. This leads to a superposition state 
\s). Recently, Thiel et al. [44] proposed a method to prepare 
all the symmetric states using a linear optical tools. We seek 
to investigate the dynamics of entanglement of the two qubits 
as they interact with the environment and with each other via 
their electric -dipoles. 

In terms of the timed Dicke basis the initial density matrix 
has only one nonzero element namely p ss (0) = 1; all the rest 
of matrix elements are zero. Since there is no initial two pho- 
ton coherence p eg (0) = 0, according to Eq. (II ldb . it remains 
zero all the time. As a consequence the expression given by 
( l23T l will be negative and hence cannot be used as entangle- 
ment measure. Moreover, it is easy to see that for initial con- 
dition we considered, p ee (t) = 0. In view of this expression 
i takes the form 



leads to 



Ci (t) = V '(fiss - Paa) 2 + 4[Im( (0as )] 2 > (24) 

and thus the concurrence can be written as 

C(i) =max(0,Ci(i)). (25) 

This expression shows that the concurrence is unity at t = 
as it should be. 

Disregarding the relative phase shift (ip = 0) the solutions 
of the elements of the density matrix in Eq. d24"l > turn out to be 

p ss (t) = exp[-2(7 + 7i 2 )4 Paa{t) = 0, and p as = 0, which 



C(t) =max(0,e- 2 ^ +712 )'). 



(26) 



We immediately see that the concurrence depends only on the 
symmetric state population, p ss (t) . As there is no single pho- 
ton coherence generated in this case, population transfer be- 
tween the \s) and \a) does not occur. As a result the initial 
entanglement experiences an enhanced decay due to the col- 
lective decay rate (712) and goes asymptotically to zero as 
t 00. For nonidentical atoms, however, even though the 
entanglement has the same behavior as identical atoms at the 
initial time, it exhibits revival at later times 111 811 . Here the de- 
tuning plays an important role in creating coherence between 
the symmetric and antisymmetric states, which is the basis for 
entanglement in the two-qubit system. In the following we 
rather show, by taking into account the spatial phase depen- 
dence of the atomic states, that quantum interference in the 
system leads to a population transfer between the symmetric 
and antisymmetric states and hence generation of coherence 

(Pas)- 

The phase shift that an atom experiences during the exci- 
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FIG. 3: (Color online) Plots of the time evolution of concurrence 
C\(t) with initial condition p ss (0) = 1, for interatomic distance 
ri2 = A/8 (712/7 = 0.9, ^12/7 — —0.9) for different values of 
£-the angle between the direction of propagation of the laser and the 
line joining the two atoms. The inset shows the populations in the 
symmetric and antisymmetric states for £ = 0(<^> = 7r/4). 
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FIG. 4: (Color online) Plots the imaginary part of coherence be- 
tween the symmetric and antisymmetric states p ils with initial con- 
dition /5 SS (0) = 1, for interatomic distance ria = A/8 (712/7 
(i.'i. Q 12 /7 = -0.9) and for £ = 0. 



tation process contain physical information about the excited 
atom. For example, in the phase factor associated with an ex- 



cited atom exp(ik • r^) the term k 



(jj h -Tj/c= L0 ti 



indicates that the atom located at position rj is excited at dif- 
ferent times ti. This has been discussed in the context of di- 
rected spontaneous emission and collective Lamb shift in re- 
cent years l29H30ll . Here we present how this phase factor can 
be used to improve the entanglement at later times. 

In one photon subspace [p ss (0) = 1] and for nonzero spa- 
tial excitation phase the important equations read l38tl 

p ss = - 2(7 + 712 cos ip)p ss - i sin 97(712 + iVLi 2 )p as 

+ i sin 93(712 -i£lu)Psa, (27) 

Paa = ~ 2(7 - 7x2 COS lfi)p aa ~ % SHI 93(712 - i£l 12 )p as 

+ i sin 99(712 + iVLi 2 )p sai (28) 

p as = - 2(7 - i£li 2 cos99)p as + i sin 93(712 + iVL 12 )p ss 

+ i sin 93(712 - i£li2)Paa- (29) 

These equations fully describe the dynamical behavior of the 
population transfer between the symmetric and antisymmet- 
ric states and the coherence developed between them under 
the given initial condition. We particularly note that the co- 
herence critically depends on the spatial excitation phase (99). 
It is not difficult to see from these equations that for a laser 
propagating perpendicular to the interqubit axis (ip — 0) there 
will be no coherence, which in turn implies the initial popula- 
tion in the symmetric state directly decays to the ground state 
without ever being transferred to the antisymmetric state. In 
this decay process the maximum entanglement present at the 
initial time will be washed out in short time. Therefore, for 
this particular initial condition, one has to play around with 
the spatial excitation phase to avoid enhanced decay of the 
entanglement. 

Using the analytical solutions of the Eqs. d27l)-(|2"9l) the con- 
currence can be expressed as 



Ci (t) =e~ 2lt [(cos 93 cosh 27^ - sinh 27i 2 t) 2 
+ sin 2 93 cos 2 2^i 2 i] 1/2 . 



(31) 



Inspection of OTI ) shows that the presence of the excitation 
phase brings in the dipole-dipole interaction (^12) into the 
dynamics. This is in contrast with the case where 93 = in 
which the concurrence is independent of the interatomic inter- 
action. Note that it is the initial preparation of the state that 
determines the dynamical behavior of the two-qubit system. 
To better understand to what extent the excitation phase mod- 
ifies the concurrence and hence the entanglement between the 
qubits, we graphically present the concurrence in Fig. [3] 

In Fig. [3] we show the evolution of the concurrence as a 
function of the angle between the direction of propagation of 
the laser and the line joining the two atoms (£) for the two- 
qubit system prepared initially in the symmetric state |s) and 
for interatomic distance ri2 = A/8. As we have discussed ear- 
lier the concurrence corresponding to ip = exhibits a sharp 
decrease and ultimately goes to zero for t — >• 00. The situa- 
tion for nonzero excitation phase is different; the concurrence 
sharply diminishes during the decay time of the symmetric 
state [27 + 712 cos 99] _1 and shows a bit of revival and decays 
slowly before it goes to zero at t — > 00. This can be under- 
stood by looking at the inset of the Fig. [3] where we plotted 
the time evolution of populations in the symmetric and anti- 
symmetric states. As can be clearly seen from this inset that 
for 99 0(£ = 0), quantum interference leads to coherent 
transfer of population from the initially populated state \s) to 
antisymmetric state \a) lf37tl and hence generation of coher- 
ence between these levels as illustrated in Fig. [4] This co- 
herence is responsible for the entanglement observed between 
the qubits. 



B. Initial mixed state 

We next consider the two qubits initially prepared in a 
mixed entangled state l24ll given by the density matrix 

p(0) = ha\l){l\ + (1 - a)|4)(4| + (b + c)|$)($|) (32) 



in which 1$) = 



:(Vb\2) 



, , „ x \fc\3)) and the normaliza- 

Jb+c v 1 ' vi// 

tion condition reads (1 + b + c)/3 = 1. Here a, b, c and x 
are independent parameters which determine the initial state 
of the two entangled qubits. Note that the above state is a 
form of generalized Werner state. The initial condition given 
by ( 132b can be written in the basis of (fT3l as 



p(0) = 



/ a 

b z 

z* c 

V 1 - a 



(33) 



C{t) = max(0,Ci(t)), 



(30) 



where z = ybc e lx is some initial single photon coherence 
in the system and \ is tne respective phase of the coherence. 
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FIG. 5: (Color online)Plots of the time evolution of concurrence C(t) 
with initial condition 6 = c = \z\ = 1 and for a = 0.6, ri2 = A/8, 
and for different values of the initial phase x 



FIG. 6: (Color online)Plots of the time evolution of concurrence with 
initial condition b = c = \z\ = 1 and for \ — t/2, ^12 = A/8 
(712/7 = 0.9, fii.2/7 = —0.9), for different values of the initial 
populations a. 



Now applying the transformation given by (l20l >. the initial 
density matrix elements for 6 = c = \z\ = 1 become 

p ee (0) = a/3,p aa (0) = (1 - cosx)/3,/9 9S (0) = (1 - o)/3 
p ss (0) = (1 + cosx)/3,p as (0) = -sinx- 

Since p es {0) = p ea (0) = /0 gs (0) = p 9a (0) = 0, the form 
of the initial density matrix remain the same, i.e., all the zero 
elements remain zero and the all the rest evolves in time. Un- 
der this scenario the expression given by d23l will be negative 
and hence cannot be an entanglement measure for the qubit 
system. Therefore, d22l is the only candidate left to quantify 
the entanglement between the qubits. For ip — the system 
of equations governing the dynamics of the two qubits can be 
solved analytically. Using these solutions, solved under the 
initial condition (|33l , the expression that describes the entan- 
glement between the qubits, C\(t), turns out to be 

Ci(t) = -e" 27 *{[(cosxcosh27i 2 t-sinh27i 2 i + a77i(t)) 2 
,1/2 



where 



sin 2 x cosh 2tt 12 i] ' - y/3a(^-m(t))}, (34) 



(2 , 2 \ 



7f 2 - 7 
27712 
7i 2 2 - 7 2 



(e- 2 ^ - cosri2 7 i 2 £), (35) 



-Aft 



2lt [ cosh 27i 2 i — cos x sinh 27i 2 £ 



+ a ( ^ +7 ' 2) (e- 2 ^-cosh2 7l2 t) 
7i2 - 7 2 

27712 



2 9 

7i 2 " 7 2 



sinh27i 2 t)] . 



(36) 



We immediately see from this result that the concurrence de- 
pends on the parameters a, which characterizes the initial pop- 
ulations of the doubly excited state i.e. when both the qubits 
are excited and on the phase parameter x which determines 
the initial populations in the symmetric and antisymmetric 
states as well as the coherence between them. If we con- 
sider that the qubits are coupled independently to their respec- 
tive vacuum environment (71 2 = 0) and are well separated 
in position (ri 2 A) so that the dipole-dipole interaction 
(f2i2 — ► 0), Ci(t) reduces to 



d(t) 



1 - sja{l - a + 2a 2 + a 4 a) 



where a(t) = \/l — exp(—2jt). Note thatCi(i) is indepen- 
dent of the initial phase y. This coincides with the earlier 
results of Yu and Eberly B2411 . 

In the following we study the dependence of the concur- 
rence and hence the entanglement between the qubits on vari- 
ous system parameters. Figure [5] shows the time evolution of 
the concurrence for ri 2 = A/8 and a = 0.6 and for different 
values of the initial phase, \. We observe from this figure that 
the initial entanglement between the qubits vanishes and ex- 
hibits revival. The amplitude of revival and the revival time 
(the time at which the entanglement revive in the system) are 
directly related to the initial coherence in the system. The 
higher the initial coherence the higher the amplitude of re- 
vival and the shorter the revival time is. Not surprisingly the 
magnitude of revival diminishes when the initial coherence 
decreases. Now keeping the initial coherence at its maximum 
value (x = tt/2), we investigate the influence of the popu- 
lation distribution between the excited and ground states on 
the concurrence. Figure [6] shows the evolution of the concur- 
rence for x = 7r/2 an d f° r different values of a. This figure 
indicates that when the initial population in the excited state 
grows the transient entanglement falls sharply and even disap- 
pears for a = 0.8 (p ee (0) ~ 0.27) in the short time window. 
The entanglement then shows revival and a slowly damping 
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the laser propagation direction and the line joining the two 
atoms. These plots clearly show that the excitation phase ef- 
fectively protects the initial entanglement from experiencing 
a sudden death and even enhances the entanglement from its 
initial value during the revival period. The amount of entan- 
glement then drops gradually and approaches zero at t oo. 
It is worth noting that the spatial excitation phase creates ad- 
ditional coherence and hence improves the revival magnitude 
over that observed for the case ip = 0. This enhanced coher- 
ence, as shown in Fig. [8j is a signature of stronger entangle- 
ment between the qubits 



yt 



V. CONCLUSION 



FIG. 7: (Color online)Plots of the evolution of concurrence with ini- 
tial condition b = c = \z\ = 1 and for ri2 = A/8 (712/7 = 
0.9, ^12/7 = —0.9), x = it/2 an d f° r the direction of propagation 
of the laser field perpendicular (£ = it/ 2) and parallel (£ = 0) to the 
line joining the two atoms. 
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FIG. 8: (Color online)Plots of the imaginary part of the one photon 
coherence with initial condition b = c = \z\ = 1 and for ri2 
A * (712/7 = 0.9, f2ia/7 = -0.9), \ = "7 2 and for the direction 
of propagation of the laser field perpendicular (£ = 7r/2) and parallel 
(£ = 0) to the line joining the two atoms. 

behavior afterwards for all values of initial populations. 

We next analyze the evolution of entanglement in the sys- 
tem by introducing the spatial excitation phase ip into the dy- 
namics. By comparing the previous results for <p = with 
the numerical plots for <p 7^ 0, we discuss the effect of the 
excitation phase on the entanglement dynamics. Our results 
are summarized in Figs. [7] and [8] In Fig. |7J we present 
a comparison of concurrence taking into account the spatial 
excitation phase ip = ir/4 (£ = 0) and in the absence of 
excitation phase, p = (£ = ir/2) for interatomic dis- 
tance less than the radiation wavelength, r 12 = A/8. Re- 
call that <p = (27r/A)ri2 cos£, where £ is the angle between 



We have investigated the effect of quantum interference in- 
duced by position dependent excitation phase on the dynam- 
ical behavior of entanglement between two interacting qubits 
coupled to a common vacuum reservoir. We have considered 
both pure and mixed entangled states for our analysis. Our 
results show that for the atoms initially prepared in a sym- 
metric state, the excitation phase induces quantum interfer- 
ence in the two-qubit system that leads to coherent popula- 
tion transfer between the symmetric and antisymmetric states. 
This thus creates a coherence which in effect slows down the 
otherwise fast decay of two-qubit entanglement considerably. 
Hence we find that the evolution of entanglement crucially 
depends on the coherence between the symmetric and anti- 
symmetric states. Furthermore, when the qubits are prepared 
in a Werner type mixed entangled state the entanglement is 
known to suffer sudden death. However, if one takes into ac- 
count the excitation phase into the dynamics the entanglement 
exhibits revival. This revival is attributed to strong coherence 
developed between the symmetric and antisymmetric states. 
A viable candidate for realization of our findings would be 
semiconductor quantum dots. Note that coupled quantum dots 
with interdot distance less than the radiation wavelength has 
already been investigated in context to quantum gates 114511 and 
photoluminescence spectra 13311 . As a future perspective, one 
can further study the effect of virtual processes on the dynam- 
ics of the system as these processes are known to influence the 
evolution of the symmetric and antisymmetric states. 
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